Surface states and the charge of a dust particle in a plasma 
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We investigate electron and ion surface states of a negatively charged dust particle in a gas 
discharge and identify the charge of the particle with the electron surface density bound in the 
polarization-induced short-range part of the particle potential. On that scale, ions do not affect 
the charge. They are trapped in the shallow states of the Coulomb tail of the potential and act 
only as screening charges. Using orbital-motion limited electron charging fluxes and the particle 
temperature as an adjustable parameter, we obtain excellent agreement with experimental data. 

PACS numbers: 52.27.Lw, 52.40.Hf, 73.20.-r, 68.43.-h 
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Motivation. — The calculation of the charge of a macro- 
scopic object in an ionized gas is one of the most funda- 
mental problems of plasma physics. It occurs in space- 
bound plasmas, where the charge of satellites is of inter- 
est jl|, in astrophysical plasmas, where one wants to know 
the charge of interstellar grains , and in laboratory gas 
discharges, where charged dust particles are either con- 
taminants, which need to be controlled, or constituents, 
whose collective properties are the subject of study [EH]- 
For laboratory plasmas, the particle charge has been 
measured in a number of experiments @, [a, 0i H 
Throughout it is thereby assumed that the particle's sur- 
face potential, and hence its charge, is the one which 
balances at the grain surface the total electron with the 
total ion charging flux: je — ji. In almost all cases, how- 
ever, the charges obtained from this condition, which is 
equivalent to forcing the net charge of the particle to 
be quasi-stationary, are too high. Usually, the approx- 
imations for the fluxes, mostly orbital motion limited 
(OML) 0, 

lull blamed for the disagreement 

and various modifications of the OML theory have been 
proposed. Although this leads sometimes to reasonable 
quantitative results 0, [^, we suspect on fundamental 
grounds that irrespective of the fluxes je = ji is not the 
condition which fixes the charge (or potential) of the par- 
ticle. 

The condition je = ji is part of the Boltzmann-Poisson 
description of the plasma-particle interaction. Its natural 
length scale is thus the length on which the Coulomb po- 
tential varies. There are however microscopic processes 
near the surface of the particle, most notable sticking 
and desorption of electrons, which affect the charge but 
take place on a much shorter scale. Once these processes 
are incorporated it is clear that the charge of the par- 
ticle is not determined by the quasi-stationarity of the 
net charge but by the individual quasi-stationarity of the 
electron and ion densities bound to the particle. This 
condition implies the former but not vice versa. It is 
thus more restrictive and leads to lower charges. 

In this Letter we describe a surface model, which ac- 



counts for plasma- and surface-induced processes, and 
calculate the charge of the particle, and its partial screen- 
ing due to trapped ions, without relying on the condition 
je ~ ji- Instead, we force the electron and ion densities 
bound to the particle to be quasi-stationary by balanc- 
ing, individually on effective surfaces, electron and ion 
charging fluxes with electron and ion desorption fluxes. 

Surface states. — We start with an investigation of 
the bound states in the static interaction potential of 
an electron (ion) and a dust particle with radius R, di- 
electric constant e, and charge —eZp. The potential 
contains a polarization-induced part, arising from the 
electric boundary condition at the grain surface, and a 
Coulomb tail due to the particle's charge [l3|- Defin- 
ing ^ = {e — l)/2{e + l)Zp, measuring distances from 
the grain surface in units of R and energies in units of 
U = ZpC^ /R, the interaction energy ai x — r/R~ 1 > x^, 
where Xb is a cut-off below which the grain surface is not 
perfect anymore, reads 
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The second line approximates the relevant parts of the 
interaction energy very well and permits an analytical 
calculation of surface states. A gas discharge usually con- 
tains enough electrons which can overcome the particle's 
Coulomb barrier U ~ o{eV). These are the electrons 
which may get bound in the polarization-induced short- 
range part of the potential, well described by the approx- 
imate expression. Ions, on the other hand, having a finite 
radius rf^^^ /R = xf^'^ « lO""*, cannot explore the poten- 
tial at these distances. The long-range Coulomb tail is 
most relevant to them, which is again well described by 
the approximate expression. 

In order to determine bound states from the 
Schrodinger equations corresponding to Ve^i(x) we have 
to specify boundary conditions. Clearly, the wavefunc- 
tions Ue,i{x) have to vanish for x ^ oo, irrespective of 
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FIG. 1: Left panel: Potential energy for an electron (ion) in 
the field of a MF particle {R = 4.7 fj,m, Z = 6800) [ll and 
s-type probability distributions shifted to the binding energy 
and maxima normalized to one. Dashed lines are the poten- 
tials used in the Schrodinger equations and thin lines are the 
bulk energy distribution functions for the hosting discharge. 
Right panel: Illustration of the surface model to be discussed 
in the main text. At quasi-stationarity, surface charges ae,i 
bound at re ~ R and « (27r(Jcxng)~^, respectively, balance 
charging fluxes Se,ije,i with desorption fluxes T~lae,i- 



the potentials. The boundary condition at Xb, in con- 
trast, depends on the potential for x < Xb, that is, on 
the surface barrier. For our purpose, it is sufficient to 
take the simplest barrier model: Ve^i{x < Xb) = oo with 



Xb 



for electrons and Xb 



for ions. The wave 



functions vanish then also at Xb and the surface states 
are basically rescaled hydrogen- type wavefunctions. 

The left panel of Fig. [1] shows s-type electron and 
ion probability densities |ue,i(a;)P (our reasoning does 
not depend on the angular momentum) for a melamine- 
formaldehyde (MF) particle {R — 4.7 /im, e = 8, and 
Zp = 6800) in a helium discharge with plasma density 
Hf, ^ rii ~ 0.62 X 10^ cm^^, ion temperature keTi = 
0.04 eV, and electron temperature ksTg = 2.2 eV 
The Rydberg series of electron surface states is only a 
few Angstroms away from the grain boundary. At these 
distances, the spatial variation of T4(a;) is comparable to 
the de-Broglie wavelength of electrons approaching the 
particle: Xf^/R IK/K'I ^ 10""*. Hence, the trapping 
of electrons at the surface of the particle is a quantum- 
mechanical effect not included in the classical description 
of the plasma-particle interaction. For ions, on the other 
hand, the lowest surface states, which carry quantum- 
mechanical features, are unimportant. Being cold and 
heavy, ions will be bound in a continuum of states be- 
low the ion ionization threshold which consists essen- 
tially of classical trapped orbits, as can be seen from the 
|Mi(a:)p for the fc (20000) state. That ions behave clas- 
sically is not unexpected. Their de-Broglie wavelength 
is much smaller then the scale on which the potential 



varies: Xf^/R - lO'^ < \V^/V{\ - 1 for x > lO'^. 

Model. — We now use the properties of the surface 
states to construct a model for the charge of the par- 
ticle. Within the sheath of the particle, the density of 
free electrons (ions) is much smaller than the density 
of bound electrons (ions). In that region, the quasi- 
stationary charge (in units of — e) is thus approximately 



Z{x) = AttR^ / dx'{l + x') 



nl{x')^n\{x') 



(2) 



with X < Xf — a/ kTi I eui , the ion Debye length, which 
we take as an upper cut-off, and ^ the density of bound 
electrons and ions. The results presented above suggest 
to express the density of bound electrons by an electron 
surface density: rig(x) w aeS{x — Xe)/ R with Xe ^ Xb ^ 
and tTg the quasi-stationary solution of 



dae/dt = Seje 



(3) 



where je is the electron charging flux from the plasma 
and Se and Tg are, respectively, the electron sticking co- 
efflcient and electron desorption time due to inelastic 
collisions between electrons and the particle (see right 
panel in Fig. [1]). We will argue below that once the par- 
ticle has collected some negative charge, not necessarily 
the quasi-stationary one, there is a critical ion orbit at 
~ 1 — 10 3> which prevents ions from hitting the 
particle surface. Thus, the particle charge is simply 



Zp = Z{Xe <X <Xi) ^ 47ri?^(sr)eje 



(4) 



For an electron to get stuck at (to desorb from) a sur- 
face it has to loose (gain) energy at (from) the surface. 
Since electrons with rather low and rather high energies 
are, respectively, reflected by the Coulomb and surface 
barrier of the particle, sticking (desorption) primarily af- 
fects electrons at energies slightly above U. Assuming 
this group of electrons to be in quasi-equilibrium with the 
surface electrons after overcoming the Coulomb barrier, 
absolute reaction rate theory [IS] allows us to estimate 
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where h is Planck's constant, Tp is the particle tempera- 
ture, and Ef is the negative of the binding energy of the 
surface state from which desorption most likely occurs. 
This phenomenological equation relates a combination of 
kinetic coefficients, which individually depend on the dy- 
namic interaction, to an energy which can be deduced 
from the static interaction alone. To go beyond Eq. ([5]) 
necessitates a quantum-kinetic treatment of the inelastic 
electron-particle interaction. 

Equation is a self-consistency equation for Zp. 
More explicitly, combined with Eq. ([5]), and using the 
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FIG. 2: Pressure dependence of the charge of a MF particle 
with R = 1 jim in the neon discharge of Ref. Q (squares). 
Solid lines denote the (isothermal) charges deduced from the 
surface model whereas dotted and dashed lines are the charges 
obtained from j^^^ = + jf'^. The inset shows the 

radius dependence of the charge for p = 50 Pa. 




FIG. 3: Radius dependence of the charge of a MF particle in 
the sheath of an argon discharge at p = 6.67 Pa \7\ (squares). 
Solid (dashed) lines give the charges deduced from the surface 
model when the depletion of Ue in the confining sheath is (is 
not) included in the OML electron charging flux. 



OML electron charging flux, which is a reasonable ap- 
proximation because, on the plasma scale, electrons are 
repelled from the particle, it reads 

Z, ^ 4.i?2_^e^e (6) 



with j^^^ = ney^k^Tj2^f^ejip[~Zpe^/RkBTe]. 
Thus, Zp depends on the radius R, the plasma parame- 
ters rie and T^, and the surface parameters Tp and E^. 

Results. — To estimate we imagine an electron with 
energy just above U approaching the grain. By neces- 
sity, it comes very close to the surface (see left panel in 
Fig. [T|). For any realistic surface barrier, the wavefunc- 
tion will therefore leak into the grain and the electron 
will strongly couple to the excitations of the grain which 
provide the thermal reservoir encoded in Tp. Hence, the 
electron will quickly relax to the lowest surface state. The 
n = 1 state for the infinitely high barrier is an approxi- 
mation to that state. Thus, E^ i?o(e - l)Vl6(e + l)^ 
where Rq is the Rydberg energy. Tp cannot be deter- 
mined as simply. It depends on the heating and cooling 
fluxes to-and-fro the grain and thus on additional surface 
parameters [31 . We use Tp therefore as an adjustable 
parameter. To reproduce, for instance, with Eq. ([6]) the 
charge of the particle in Fig. [1] Tp = 395 K implying 
E^ « 0.51 eV and (st)^ « 0.4 x IQ-^ g. 

In Fig. [5] we analyze within our approach the pres- 
sure dependence of the charge of a MF particle with 
R ~ 1 fim in the bulk of the neon discharge of 
Ref. Since the plasma parameters entering Eq. ^ 
are known 0, ?p is again the only free parameter. Fix- 
ing Tp at a particular value gives the isothermal particle 
charges Zp(Tp). From Zp{Tp) = Z^xp follows then the 



Tp required to reproduce the data. The predicted in- 
crease of Tp with pressure is realistic. Indeed, assum- 
ing Tp ~ Tg, with Tg the gas temperature, Tp is in 
accordance with what one would expect from the pres- 
sure dependence of Tg in noble gases For com- 
parison we also plot the particle charges deduced from 

^■OML ^ ^-OML ^ ^-CX ^iti^ ^-OML ^ n,^ ksT, /2nm,[l + 

Zpe^/RksTi] the OML ion charging flux and jf^ = 
n^iO.lXf /lcK)^/kBT^/2^^m^{Zpe^ /RksT)^ the ion flux 
due to charge-exchange (CX), where /cx — (f^cxng)^^ 
is the scattering length, ctcx is the cross section, and 
rig — p/kBTg is the gas density Although for = 
2 X lO^^** cm^ (which may be unrealistically large [17|) 
the agreement with the data is good, the radius depen- 
dence of Zp at fixed pressure shown in the inset indicates 
that something must be wrong with the flux balance cri- 
terion. The data clearly appear to be closer to the non- 
linear i?— dependence obtained from the surface model 
than to the linear one resulting tromj^^^ — ji*^^ + ji^'^ ■ 
Figure [3l showing the i?— dependence of Zp for MF 
particles conflned in the sheath of an argon discharge 0] , 
provides additional support for our model. To approxi- 
mately account for the fact that particles with different 
radius experience different plasma environments, we in- 
cluded the depletion of rie in the sheath by replacing rie in 
^■OML neexp[e$(zeq(i?))/fcsTe] with $(z) the sheath 
potential and Zeq(-R) the equilibrium position of the par- 
ticle with radius R [7|] . When the grains are not too deep 
in the sheath (i? < 5 fim), we flnd excellent agreement 
with the data for Tp = 420 K. Our approach fails, how- 
ever, for i? > 5 /im (see inset). We attribute this to the 
ad-hoc description of je which may not capture the total 
electron charging flux close to the electrode. 
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Equation ([6]) depends on the assumption that once the 
particle is negatively charged ions are trapped far away 
from the grain surface. Indeed, a recent study based on 
the Boltzmann-Poisson equations has shown that charge- 
exchange collisions lead to a local pile-up of ions in the 
sheath of the particle [3, We come to the same 

conclusion from the surface physics point of view. Sim- 
ilar to an electron, an ion gets bound to the grain only 
when it looses energy. Because of the long-range at- 
tractive ion-grain interaction, the ion will be initially 
bound far away from the grain surface (see left panel 
in Fig. [T]) . The coupling to the excitations of the grain is 
thus negligible and only inelastic processes due to the 
plasma are able to induce transitions to lower bound 
states. Since the interaction is classical, inelastic col- 
lisions, for instance, charge-exchange between ions and 
atoms, act like a random force. Energy relaxation can be 
thus envisaged as a de-stabilization of orbits whose spa- 
tial extension is comparable to or larger then the scat- 
tering length. Smaller orbits are unaffected because the 
collision probability during one revolution is vanishingly 
small. For a circular orbit, a rough estimate for the criti- 
cal radius is = R{l + Xi) = (27r(Tcx'T-g)~^ which leads to 
Xi ^ 5 ^ ^ when we use the parameters of the he- 
lium discharge of Fig. [T] and CTcx = 0.32 x IQ-^^ cm^ 
Thus, there is a relaxation bottleneck at Xi and ions are 
trapped deep in the sheath of the particle. 

To determine the partial screening due to trapped ions 
we model the ion density accumulating in the vicinity 
of the critical orbit by a surface density ai which balances 
at Xi the ion charging flux with the ion desorption flux 
(see right panel in Fig. [T|). Mathematically, this gives 
rise to a rate equation similar to Q but now for the 
ions. At quasi-stationarity, the ion surface density is thus 
o'i = {sT)iji. Although Eq. ([5]) assumes excitations of the 
grain to be responsible for sticking and desorption we 
expect a similar expression (with , Tp replaced by Ef, 
Tg) to control the density of trapped ions. From Eq. ^ 
we then obtain Z{xi < x < Xf) = Zp — Zi with 



Z, = 47ri?2(l 



h 



,^Et{Zp)/kBT,jB 



(7) 



the number of trapped ions. Since the critical orbit is 
near the particle-sheath-plasma boundary, it is fed by 
the Bohm ion flux = O.Gni-^fcsTe/mi. The ion 
desorption energy is the negative of the binding en- 
ergy of the critical orbit, Ef{Zp) — -Vi{xi)U{Zp) = 
ATTacx.o.BngZpRo, and depends strongly on Zp and Xi. For 
the situation shown in Fig.[l]we obtain Ef « 0.37 eV and 
(sr), « 0.6 X 10"* s when we use Tg = Tp = 395 K, the 
particle temperature which reproduces Zp « 6800. The 
ion screening charge is then Zi « 148 <§; Zp which is the 
order of magnitude expected from molecular dynamics 
simulations 2l|. Thus, even when the particle charge is 
defined hy Z{xi < x < \f) it is basically given by Zp. 



Summary. — We constructed a surface model to cal- 
culate the charge (partial screening) of a particle in a 
plasma by balancing, on an effective surface, the elec- 
tron (ion) charging with the electron (ion) desorption 
flux. The number of electrons bound in the polariza- 
tion potential determines the charge of the particle. Us- 
ing the grain temperature as an adjustable parameter 
we obtained far better agreement with measurements, in 
particular, with respect to the radius dependence of the 
charge, then approaches based on balancing at the grain 
surface the total charging fluxes which we argue is the 
wrong condition. It neglects the microscopic processes 
determining the charge of the particle: sticking and des- 
orption of electrons at the grain surface. 
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